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Abstract 

This paper deals with the numerical computation of boundary null controls for the ID 
wave equation with a potential. The goal is to compute approximations of controls that drive 
the solution from a prescribed initial state to zero at a large enough controllability time. We 
do not apply in this work the usual duality arguments but explore instead a direct approach 
in the framework of global Carleman estimates. More precisely, we consider the control that 
minimizes over the class of admissible null controls a functional involving weighted integrals 
of the state and the control. The optimality conditions show that both the optimal control 
and the associated state are expressed in terms of a new variable, the solution of a fourth- 
order elliptic problem defined in the space-time domain. We first prove that, for some specific 
weights determined by the global Carleman inequalities for the wave equation, this problem 
is well-posed. Then, in the framework of the finite element method, we introduce a family 
of finite-dimensional approximate control problems and we prove a strong convergence result. 
Numerical experiments confirm the analysis. We complete our study with several comments. 

Keywords: one-dimensional wave equation, null controllability, finite element methods, Carleman 
estimates. 
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1 Introduction. The null controllability problem 

We are concerned in this work with the null controllability for the ID wave equation with a 
potential. The state equation is the following: 

Vtt ~ {a(x)y x ) x + b(x, t)y = 0, (x, t) e (0, 1) x (0, T) 

y(0,t)=0, y(l,t)=v(t), te(0,T) (1) 

y(x,0) = y Q (x), y t (x,0) = yi(x), x e (0,1). 

Here, T > and we assume that a E C 3 ([0, 1]) with a(x) > a > in [0, 1], b E L°°((0, 1) x (0,T)), 
yo E L 2 (0, 1) and y\ E if _1 (0, 1); v — v(t) is the control (a function in L 2 (0,T)) and y = y(x,t) is 
the associated state. 

In the sequel, for any r > we denote by Q T and S r the sets (0, 1) x (0, r) and {0, 1} x (0, r), 
respectively. We will also use the following notation: 

Ly:=y tt - (a(x)y x ) x + b(x,t)y. (2) 
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For any (yo,Vi) G Y := L 2 (0, I) x H 1 (0, 1) and any t> G L 2 (0,T), it is well known that there 
exists exactly one solution y to (JlJ, with the following regularity: 

yeC ([0,T];L 2 (0,l))nC 1 ([0,T]; J ff- 1 (0,l)) (3) 

(see for instance (2"3"]V 

On the other hand, for any T > 0, the null controllability problem for (JlJ at time T is the 
following: for each (t/o>2/i) £ Y, find u 6 L 2 (0,T) such that the corresponding solution to |l| 
satisfies 

y(-,T) = 0, y t (-,T)=0 in (0,1). (4) 

In view of the linearity and reversibility of the wave equation, (JlJ is null-controllable at T if 
and only if it is exactly controllable in Y at time T, i.e. if and only if for any (yo,yi) G Y and 
any (zq, z±) E Y there exist controls v E L 2 (0, T) such that the associated y satisfies 

y(',T)=Zo, y t (-,T) = z 1 in (0,1). 

It is well known that is null-controllable at any large time T > T* for some T* that depends 
on a (for instance, see [2] [23J for a = 1 and 6 = leading to T* = 2 and see [32] for a general 
situation). As a consequence of the Hilbert Uniqueness Method of J.-L. Lions [23] . it is also known 
that the null controllability of Q is equivalent to an observability inequality for the associated 
adjoint problem. 

The goal of this paper is to design and analyze a numerical method allowing to solve the 
previous null controllability problem. 

So far, the approximation of the minimal L 2 -norm control — the so-called HUM control — has 
focused most of the attention. The earlier contribution is due to Glowinski and Lions in [19] (see 
also [21j for an update) and relies on duality arguments. Duality allows to replace the original 
constrained minimization problem by an unconstrained and a priori easier minimization (dual) 
problem. However, as observed in [T5] and later in [31], depending on the approximation method 
that is used, this approach can lead to some numerical difficulties. 

Let us be more precise. It is easily seen that the HUM control is given by v(t) — a(l)4> x (l,t), 
where solves the backwards wave system 

L<j> = in Qt 

(f> = on E T (5) 

(<K-,T),M',T)) = (<h,<M in (0,1) 

and (0o,</>i) minimizes the strictly convex and coercive functional 

1 f 1 

^(0o, 0i) = -11^(1)0^(1, OlllaccT) + J ya{x) 4>t{x,Q) dx - (yi,<j>{- ,0))h-i,h* (6) 

over H = ^(0,1) x L 2 (0,1). Here (-,-) H -i H i denotes the duality product for i/ _1 (0, 1) and 
The coercivity of I over H is a consequence of the observability inequality 

ll0o||ffi( O ,i) + II</>i|Il2(o,i) <C\\<Px(l,-)\\ 2 L 2 {QX) Vf>o,<Ai) e H, (7) 

that holds for some constant C = C(T). This inequality has been derived in [23] using the 
multipliers method. 

At the numerical level, for standard approximation schemes (based on finite difference or fi- 
nite clement methods), the discrete version of ([7| may not hold uniformly with respect to the 
discretization parameter, say h. In other words, the constant C — C(h) may blow up as h goes 
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to zero. Consequently, in such cases the functionally (the discrete version of X) fails to be coercive 
uniformly with respect to h and the sequence {vh}h>o may not converge to v as h — > 0, but diverge 
exponentially. These pathologies, by now well-known and understood, are due to the spurious 
discrete high frequencies generated by the finite dimensional approximation; we refer to |34| for a 
review on that topic; see [23] for detailed examples of the behavior observed with finite difference 
methods. 

Several remedies based on more elaborated approximations have been proposed and analyzed 
in the last decade. Let us mention the use of mixed finite elements 6J, additional viscosity terms 
which have the effect to restore the uniform property [TJ [23] and also filtering technics [TTj . Also, 
notice that some error estimates have been obtained recently, see [7[ QT] . 

In this paper, following the recent work |12j devoted to the heat equation, we consider a different 
approach. Specifically, we consider the following extremal problem: 

Minimize J{ Vl v)= l -JJ^ p 2 \y\ 2 dx dt + l - jf pl\v\ 2 dt 
Subject to (y,v) € C(y ,yi;T) 
where C(yo, y±; T) denotes the linear manifold 

C(y ,yi\T) = { (y,v) : v e L 2 (0,T), y solves and satisfies Q }. 

Here, we assume that the weights p and po are strictly positive, continuous and uniformly 
bounded from below by a positive constant in Qj* and (0,T), respectively. 

As in the previous L 2 -norm situation (where we simply have p = and po = 1), we can apply 
duality arguments in order to find a solution to Q , by introducing the unconstrained dual problem 

Minimize J* (m,<£o,<M = J ff p~ 2 \p\ 2 dx dt +\ f p^ 2 |a(l)0 x (l, t)\ 2 dt 

1 JJQt 1 JO 

1 (9) 
ya(x)<f>t{x,0)dx- {yi, 0(- , 0)) H -i H i 







, Subject to (fj,, 0o, <Pi) € L 2 (Q T ) x H, 

where solves the nonhomogeneous backwards problem 

Lcj) = fj, in Qt 

= on 
(0(-,T),0 t (.,T)) = (0o,0i) in (0,1). 

Here, J* is the conjugate function of J in the sense of Fenchel and Rockafellar JTHl HH] and, 
if p G L°°(Qt) and pa 6 L°°(0, T) (that is, p~ 2 and p^ 2 are positively bounded from below), J* is 
coercive in L 2 (Qt) x H thanks to Therefore, if (fx, 0Oi0l) denotes the minimizer of J*, the 
corresponding optimal pair for J is given by 

v = -a(l)pQ 2 tp x (l, •) in (0,T) and y = -p~ 2 p in Q T . 

At the discrete level, at least for standard approximation schemes, we may suspect that the 
coercivity of J* may not hold uniformly with respect to the discretization parameters, leading to 
the pathologies and the lack of convergence we have just mentioned. 

On the other hand, the fact that the state variable y appears explicitly in the cost J makes it 
possible to avoid dual methods. We can use instead suitable primal methods to get an optimal 
pair (y,v) € C(yo,yi]T). The formulation, analysis and practical implementation of these primal 
methods is the main goal of this paper. 
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More precisely, the optimality conditions for the functional J allow to express explicitly the 
optimal pair (y, v) in terms of a new variable, the solution of a fourth-order elliptic problem in the 
space-time domain Qt that is well-posed under some conditions on T, the coefficient a and the 
weights p and po- Sufficient conditions are deduced from an appropriate global Carleman estimate, 
an updated version of the inequalities established in [3]. From a numerical viewpoint, this elliptic 
formulation is appropriate for a standard finite element analysis. By introducing adequate finite 
dimensional spaces, we are thus able to deduce satisfactory convergence results for the control, 
something that does not seem easy to get in the framework of a dual approach. 

A similar primal approach, based on ideas by Fursikov and Imanuvilov [16] , has been used 
in |12j for the numerical null controllability of the heat equation. 



This paper is organized as follows. 

In Section [2j adapting the arguments and results in [T^], we show that the solution to ^ can 



be expressed in terms of the unique solution p to the variational problem ( 23 1 in the Hilbert space 



P, defined as the completion of Pq with respect to the inner product ( 18 ); see Proposition 2.2 The 



well-posedness is deduced from the application of Riesz's Theorem: a suitable global Carleman 



inequality ensures the continuity of the linear form in ( 23 1 for T large enough when p and po are 



given by (Tl9|); see Theorem 2.1 



In Section [3j we analyze the variational problem ( 23 ) from the viewpoint of the finite element 
theory. Thus, we replace P by a conformal finite element space Ph of C 1 (Qt) functions defined 



by ( 33 ) and we show that the unique solution p^ £ P^ to the finite dimensional problem ( 38 ) 
converges (strongly) for the P-norm to p as h goes to zero. 

Section [4] contains some numerical experiments that illustrate and confirm the convergence of 
the sequence {ph}- 

Finally, we present some additional comments in Section [5] and we provide some details of the 



proof of Theorem 2.1 in the Appendix. 



2 A variational approach to the null controllability problem 

With the notation introduced in Section [TJ the following result holds: 

Proposition 2.1 Let T > be large enough. Let us assume that p and po are positive and satisfy 
p £ C°(Qt), Po £ C°(0,T) and p, po > p > 0. Then, for any (yo,Ui) S Y , there exists exactly one 
solution to the extremal problem 

The proof is simple. Indeed, for T > T* , null controllability holds and C(yo, yi; T) is non-empty. 
Furthermore, it is a closed convex set of L 2 (Qt) x L 2 (0,T). On the other hand, (y,v) h-> J(y,v) 
is strictly convex, proper and lower-semicontinuous in L 2 (Qt) x L 2 (0,T) and 

J(y,v) +oo as \\(v,v)\\ L 2(q t)xL 2(x t ) -> +oo. 

Hence, the extremal problem ^ certainly possesses a unique solution. 

In this paper, it will be convenient to assume that the coefficient a belongs to the family 

A{x a ,a ) = {a £ C 3 ([0, 1]) : a{x) > a >0, 

— min (a(x) + (x — xo)a x (x)) < min (a(x) + — (x — xo)a x (x) \ } ^ 
[0,1] \ / [0,1] \ 2^ / 

where Xo < and ag is a positive constant. 
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It is easy to check that the constant function a(x) = ao belongs to A(xo, ao). Similarly, any non- 
decreasing smooth function bounded from below by ao belongs to A(xq,oo). Roughly speaking, 
a G A(xq, ao) means that a is sufficiently smooth, strictly positive and not too decreasing in [0, 1]. 



Under the assumption ( 10 ), there exists "good" weight functions p and po which provide a very 
suitable solution to the original null controllability problem. They can be deduced from global 
Carleman inequalities. 

The argument is the following. First, let us introduce a constant (3, with 

— mm (a(x) + (x — Xo)a x (x)j < (3 < mm (a{x) + — (x — Xo)a x (x) S j (11) 
and let us consider the function 

</>(x,t) := \x ~ x \ 2 - pt 2 + M , (12) 

where Mq is such that 

<f>{x, t) > 1 V{x, t) e (0, 1) x (-T, T), (13) 
i.e. M Q > 1 - \x Q \ 2 + j3T 2 . Then, for any A > we set 

<p(x,t) : =e A ^*). (14) 

The Carleman estimates for the wave equation are given in the following result: 

Theorem 2.1 Let us assume that xq < ; ao > and a G ^4(a;o,ao)- Let (3 and ip be given 



respectively by (11) and (14 1. Moreover, let us assume that 

T > - max a(x) 1/2 (x - x ). (15) 
P [o,i] 

Then there exist positive constants so and M , only depending on Xq, oq, ||a||c 3 ([0,l])j ll^lli 00 (Qr) 
and T , such that, for all s > Sq, one has 

s[ [ e 2sv (\w t \ 2 + \w x \ 2 ) dx dt + s 3 [ [ e 2sip \w\ 2 dx dt 

f T ,1 f T ( 16 ) 

<M / e 2sip \Lw\ 2 dxdt + Ms e 2stp \w x (l, t)\ 2 dt 
J-tJo J-t 

for any w G L 2 (-T, T; Hq(0, 1)) satisfying Lw G i 2 ((0, 1) x (— T, T)) and w x (l,-) G L 2 {-T,T). 

There exists an important literature related to (global) Carleman estimates for the wave equa- 
tion. Almost all references deal with the particular case a = 1; we refer to [31 HI HHJ, [23 S3] • The 
case where a is non-constant is less studied; we refer to 1171. 



The proof of Theorem 2.1 follows closely the ideas used in the proofs of Theorems 2.1 and 2.5 
in [4] to obtain a global Carleman estimate for the wave equation when a = 1. The parts of the 
proof which become different for non-constant a are detailed in the Appendix of this paper. 

In the sequel, it is assumed that xq < and ao > are given, a G A(xq, ao) and 
2 



T > ^ max a(x) 1/2 (a; - x ), with f3 satisfying ([IT]). (17) 



Let us consider the linear space 

Po = {qeC co (Q^) :g = 0onE T }. 
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The bilinear form 



p 2 LpLqdxdt 



Qt 



Po 2 a(l) 2 p x {l,t)q x (l,t)dt 



(18) 



is a scalar product in Pq. Indeed, in view of (17), the unique continuation property for the wave 



equation holds. Accordingly, if q 6 P 0l Lq = in Qt and q x — on {1} x (0, T), then q = 0. This 
shows that (• , -)p is certainly a scalar product in Pq. 

Let P be the completion of P with respect to this scalar product. Then P is a Hilbert space 



for (• , -)p and we can deduce from Theorem 2.1 the following result, that indicates which are the 
appropriate weights p and po for our controllability problem: 



Lemma 2.1 Let us assume that s > so, let us set 

p{x,t) := e- s ^ 2t - T \ p (t) := p(l,t) 



(19) 



and let us consider the corresponding Hilbert space P. Then there exists a constant Co > 0, only 
depending on Xq, ao, ||<j||c 3 ([q,i])j IHIl=°(q t )' ^> s andT, such that 



lb(" >°)Hh Mo,i) + \\Pt(- ,0)\\h{o,i) ^ c ° (p>p)p e p - 

Proof: For every p € P, we denote by p e L 2 ((0, 1) x (— T, T)) the function defined by 

t + T 



(20) 



p(x,t) =p[x 



It is easy to see that p e L 2 (-T,T; H^(fl)), Lp £ £ 2 ((0,1) x (-T, T)) and p x (l,-) G L 2 (-T,T), 
so that we can apply Theorem |2.1| to p. Accordingly, we have 



si I e 

-T JO 



\p x \ 2 ) dx dt + .s- 



' /T f 1 e 2sv \p\ 2 dxdt 

T JO 



< C 



2 * v (\Pt 

f f e 2sv \Lp\ 2 dxdt + Cs f e 2s ^ x ^\p x {l,t)\ 2 dt 

J-TJO J-T 



(21) 



where C depends on x , a , ||a|| C 3([o,i]) ; II^IIl~(q t ) and T. 



following for any T satisfying (17) 



Replacing p by its definition in (21 ) and changing the variable t by t' = 2t — T we obtain the 

p- 2 {\p t \ 2 + \p x \ 2 )dxdt + s 3 [[ p- 2 \p\ 2 dxdt 



<C [[ p- 2 \L P \ 2 dxdt + Cs [ Po 2 \p x (l,t)\ 2 dt, 
JJqt Jo 

where C is replaced by a slightly different constant. Finally, from Corollary 2.8 in £3], we obtain 



the estimate ( 20 ) 



Remark 1 The estimate ( 20 ) must be viewed as an observability inequality. As expected, it holds 



if and only if T is large enough. Notice that, when a(x) = 1, the assumption (17) reads 

T>2(l-x ). 



This confirms that, in this case, whenever T > 2, (20) holds (it suffices to choose xo appropriately 

□ 



and apply Lemma 2.1 see |23j ) 
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The previous results lead to a very useful characterization of the optimal pair (y, v) for J: 



Proposition 2.2 Let us assume that s > Sq, let us set p and p Q as in (19) and let us consider the 
corresponding Hilbert space P. Let (y,v) £ C(yo,yi,T) be the solution to Q. Then there exists 
p £ P such that 

V = -p~ 2 Lp, v = -(a(x)po 2 p x )\ x=1 . (22) 
Moreover, p is the unique solution to the following variational equality: 

pT 



Qt 



p z LpLqdxdt+ I p 2 a 2 (l)p x (l, t) q x (l, t) dt 



= y o (x)q t (x,0)dx- (y 1 ,q(-,0))H- 1 ,m Vq e P; p € P. 
Jo 



(23) 



Here and in the sequel, we use the following duality pairing: 

(y\q(;0)) H -i,HZ = J o ^((-A)- 1 y 1 )(x)q x (x,0)dx, 
where — A is the Dirichlet Laplacian in (0, 1). 



Proof: From the definition of the scalar product in P, we see that p solves (23) if and only if 
(P,q)p= yo(x)q t (x,0)dx - (y 1 ,q(-,0)) H -i H i Vq £ P; p £ P. 



In view of Lemma 2.1 and Riesz's Representation Theorem, problem (23) possesses exactly one 
solution in P. 



Let us now introduce y and v according to (22) and let us check that (y,v) solves ([8]). First, 
notice that y £ L 2 (Qt) and v £ L 2 (0,T). Then, by replacing y and v in (23), we obtain the 
following: 



yLqdxdt+ / a(l)v{t)q x (l,t)dt= / y {x) q t (x, 0) dx - (y 1 , q(; 0)) H -i m Vq £ P. (24) 



Hence, (y, v) is the solution of the controlled wave system ([!]) in the transposition sense. Since 
y £ L 2 (Q T ) and v £ L 2 (0,T) the couple (y, v) belongs to C(y ,yi,T). 

It remains to check that (y,v) minimizes the cost function J in But this is easy. Indeed, 
for any (z, w) £ C(yo,yi,T) such that J(z,w) < +oo, one has: 

J(z,w) > J(y,v) + // p 2 y (z — y) dx dt + / p 2 ] v{w — v)dt 
JJqt Jo 

= J(y,v) — // Lp(z — y)dxdt+ / Pqv(w — v) dt — J(y, v). 
JJqt Jo 

The last equality follows from the fact that 



Lp(z — y)dxdt= // pL(z — y)dxdt 
<)i JJqt 

1 

T 



+ / [Pt (z - y)} dx - [< (z - y) t ,p >H-i,ffi] 



[a(x)p x (z - y)]l dt + / [a(x)p (z - y) x \l dt 



the boundary condition for p (see Remark [2] below), the fact that both (y,v) and (z,w) belong 
to C(y , yi ;T) and p}. □ 
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Remark 2 From (22) and (23), we see that the function p furnished by Proposition 2.2 solves, at 
least in the distributional sense, the following differential problem, that is of the fourth-order in 
time and space: 

' L(p- 2 Lp) = 0, {x,t)eQ T 
p(0,t) = 0, (p- 2 Lp)(0,t) =0, te(o,T) 
p(M)=0, (p-'Lp + a/^VXM) =0, te(o,r) 
(/r 2 Lp)(z,0)=2/o(z), (p 2 Lp)(x, T) — 0, as e (0, 1) 
[ (p- 2 ip) t (ar,0) = 2/ 1 (x), (p- 2 J Lp) t (a;,T)=0, ie(0,l). 



(25) 



Notice that the "boundary" conditions at t = and t = T are of the Neumann kind. 



□ 



Remark 3 The weights p 1 and p 1 behave exponentially with respect to s. For instance, we 
have 

P (x,t)- X =ex P ( S e A d— o| 2 ^(2 t -T) 2 +M ) 



For large values of the parameter s (greater than so > 0, see the statement of Theorem 2.1 ), the 
weights p~ 2 and Pq 2 may lead in practice to numerical overflow. One may overcome this situation 
by introducing a suitable change of variable. 

More precisely, let us introduce the variable z = pp and the Hilbert space M — pP, so that the 
formulation (23) becomes: 



Qt 



p 2 L(pz)L(pz)dxdt + / p 2 a 2 {l){pz) x {l,t)(pz) x (l,t)dt 
Jo 

y o (x)(pz) t {x,0)dx-( yi ,(pz)(-,0)) H -i H i Vz G M; z G M. 



(26) 



The well-posedness of this formulation is a consequence of the well-posedness of ( 23 1 . Then, after 
some computations, the following is found: 

pL(p~ 1 z) = p^ 1 (ipz) t - (a(pz) x ) x + bpz^j 

= (p~ 1 p t )z + z t - a x ((p~ 1 p x )z + z x ) - a(2p~ 1 p x z x + p^ 1 p xx z + z xx ) + bz 



with 

P~ l Px 

Similarly, 



-stp x (x,2t-T), p 1 p t = ~2sip t (x,2t -T), p 1 p xx = -sip xx + (stp x )* 
(Po V)*)(M) =z x {l,t). 



Consequently, in the bilinear part of ( 26 1 , there is no exponential (but only polynomial) function 
of s. In the right hand side (the linear part), the change of variable introduces negative exponentials 
in s. A similar trick has been used in [12 in the context of the heat equation, where we find weights 
that blow up exponentially as t — > T~ . □ 

Remark 4 The exponential form of the weights p and po is purely technical and is related to 
Carleman estimates. Actually, since for any s and A these weights are uniformly bounded and 
uniformly positive in Qt-, the space P is independent of p and po and one could apply the primal 
approach to the cost J (defined in (J8|) for any bounded and positive weights. In particular, one 
could simply take p = 1 and po = 1; the estimates (20) would then read as follows: 

lb(-.°)llHj(o,i) + WM-, 0)111.(0,1) < Co (||ip|| 2 L2(QT) + ||a(l)p :c (l,-)l| 2 L 2 (o,T)) VpeP (27) 

for some constant Cq > 0. This inequality can also be obtained directly by the multipliers method; 
we refer to E2 and references therein. □ 
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Remark 5 As remarked in [3] (see Remark 2.7), the estimate (20) can be proven for a weight 
po which blows up at t — and t = T. For this purpose, we consider a function 9g £ C 2 ([0,T]) 
with 9$(0) — 9g(l) = and 9$(x) — 1 for every x G (5,T — S). Then, introducing again p(x,t) := 



9s(t)p(x, (t + T)/2), it is not difficult to see that the proofs of Lemma 2.1 and Theorem 2.1 can be 



adapted to obtain ( 20 ) with 

p(x,t) = e - s ^ 2t - T \ 



Po(t) = 5 (t)- 1 p(l,t). 



Thanks to the properties of 9s, the control v defined by 

vanishes at t = and also at t = T, a property which is very natural and useful in the boundary 
controllability context. In the sequel, we will use this modified weight po, imposing in addition, 
for numerical purposes, the following behavior near t = and t = T: 



lira m = 0(1), lim -P= = C9(1). 



t-)-0+ y/t 



(28) 

□ 



3 Numerical analysis of the variational approach 



We now highlight that the variational formulation ( 23 1 allows to obtain a sequence of approxima- 



tions {vh} that converge strongly towards the null control v furnished by the solution to ([8]) 

3.1 A conformal finite dimensional approximation 

Let us introduce the bilinear form m(-, •) over P x P 

m(p,q) := (p,q) P = // p~ 2 LpLqdxdt+ / a(l) 2 p^ 2 p x (l,t) q x (l,t) dt 
JJqt Jo 

and the linear form £, with 

{£,q) ~ / yo(x)q t {x,0)dx - (y 1 ,q(-,0)) H -i H i. 



Then (23 1 reads as follows: 



m(p,q) = (e,q), VgeP; p G P. 



(29) 



Let us assume that a finite dimensional space Ph C P is given for each h £ K^_. Then we can 
introduce the following approximated problems: 



m(ph,qh) = (i,qh), Vgy, G P h ; Ph G Ph- 



(30) 



Obviously, each (30) is well-posed. Furthermore, we have the following classical result: 



Lemma 3.1 Let p G P be the unique solution to (29) and let p^ £ P^ be the unique solution to 



(30). Then we have: 



\\p-Ph\\p < in t \\p-qh\\p- 
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Proof: We write that 

\\Ph -p\\p = m(Ph ~P,Ph ~P) = m{Ph ~P,Ph - Qh) +m(p h -p,q h -p)- 

The first term vanishes for all qh £ Ph- The second one is bounded by \\ph — p\\p\\qh ~ p\\p- So, 
we get 

\\p-Ph\\p <\\p-qh\\p Vqh e Ph 

and the result follows. □ 

As usual, this result can be used to prove that ph converges towards p when the spaces Ph are 
chosen appropriately. More precisely, let us assume that an interpolation operator II ^ : Pq — > Ph 
is given for any h 6 and let us suppose that 

\\p-U h p\\p^0 as h-> (0,0) VpeP . (31) 

We then have the following convergence result: 



Proposition 3.1 Let p £ P be the solution to (29) and let p h £ Ph be the solution to ( 30 ) for 
each h £ R+. Then 

\\p-Ph\\p^0 as h^ (0,0). (32) 

Proof: Let us choose e > 0. Since Po is dense in P, there exists p e £ Pq such that \\p — Pe\\p < e - 
Therefore, we find from Lemma HO] that 



\\P-Ph\\p < \\p - UhPe\\p 

< \\p-Pe\\p + \\Pe - TlhPe\\p 

<e + \\p e - U h pe\\p- 



But we know from (31) that \\p e — n^pe||p goes to zero as h £ Ki, h — > (0,0). Consequently, we 



also have ( 32 ) . □ 



3.2 The finite dimensional spaces Ph 

The spaces Ph must be chosen such that p~ x Lph belongs to L 2 (Qt) for any Ph £ Ph- This means 
that ph must possess second-order derivatives in L^ oc (Qt)- Therefore, a conformal approximation 
based on a standard quadrangulation of Qt requires spaces of functions continuously differentiable 
with respect to both variables x and t. 

For large integers N x and N t , we set Ax = 1/N X) At = T/N t and h — (Ax, At). We introduce 
the associated quadrangulations Qh, with Qt = V)k^q h ^ an< ^ we assume that {Qh}h>o is a 
regular family. Then, we introduce the space Ph as follows: 

Ph = { z h £ C\Qt) : z h \ K £ F(K) £ Q h , z h =0 on S T }. (33) 

Here, ¥(K) denotes the following space of polynomial functions in x and t: 

F(K) = (F 3 , x ®F 3 , t )(K) (34) 

where P r ^ is by definition the space of polynomial functions of order r in the variable £. 
Obviously, Ph is a finite dimensional subspace of P. 
Let us introduce the notation 



K k i := [x fc ,Xfc+i] x [ti,ti+i], 
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where 

x k := (k - 1) Ax, ti := (I - 1) At, for k = 1, . . . , N x + 1, I = 1, . . . , N t + 1. 

For any fc, we denote by (-£<ifc)o<i<3 the Hermite functions associated to [x kl x k+ i]. They are 
given by 

Lofc(z) := (1 + 2c)(l - c) 2 , L lfe (x) := c 2 (3 - 2c) 
£-2fc(z) := Axc(l - c) 2 , L 3 fe(^) := Axc 2 (c - 1) 
c := (x - x k )/Ax. 
Recall that, for any / e C 1 ([x k , Xfc+i]), the function 

l l 
(n A x/)(a;) := ^iifc(a;)/(x i+fc ) + L l+2 , k (x)f x (x l+k ) 

i=0 i=0 

is the unique element in P 3 ([xfe, x/c+i]) that satisfies 

(n A x/)(^fe+i) = f{Xk+i), (ttAxf)x(xk+i) = {fx)(xk+i), i = 0, 1. 

In a similar way, we denote by {£ji)o<j<3 the Hermite functions associated to the time in- 
terval [ti,ti + i\. Then, from the definition of ¥(K k[ ), we can obtain easily for any u £ P the 
polynomial function in F(K k i) uniquely determined by the values of u, u x , u t and u xt at the 
vertices of K k i : 

Lemma 3.2 For each u e P , let us define the function TLhU as follows: for any k and I, 
l l 
n h u(x,t) := ^2 L ik(x)£ji(t)u(x i+k ,t j+l ) + ^2 L i+2,k(x)C 1 i(t)u x (x l+k ,t J+ i) 

i,j—0 i>j=Q 

1 1 

+ ^2 L ik {x)C j+ 2d{t)u t {x l+k ,t j+ i) + ^2 L i+2 . k {x)C : j + 2j{t)u xt {x i+k ,t j+ i) 



in K H = [x k ,x k + Ax] x [t u ti + At]. 

Then H h u is the unique function in Ph that satisfies 

U h u(x k+i ,tl + j) = u(Xk+i,tl + j), (U h u(Xk+i,tl + j)) x = U x {Xk+i,tl+j), 
(Il h u(Xk + i,tl + j)) t = Ut(Xk+i,tl+j), (U h u(Xk+i,tl + j)) xt = U xt (x k +i,tl + j) 

for all i,j G {0, 1}. The linear mapping YLh : P h-» Ph is by definition the interpolation operator 
associated to Ph- 

This result allows to get an expression of u — UhU on each element K k i that will be used in the 
next section: 

Lemma 3.3 For any u <G Pq, we have 

u- U h u = ^2 [mijU x (x i+k ,t j+ i) +nijU t (x i+ h,tj + i) + Pij u tx (x i+k ,t j+l )) 

i i'=0 ^ ' 

(35) 

+ ^ L ik £jiU[u;x l+k ,tj + i} 

i,3=0 

in K k i , where the , and Pij are given by 

mij(x,t) := (^L lk (x)(x - - L i+2 ,k(x)jjCj(t) 

nij(x,t) := L ik (x) ^Cj(i)(t - t 3 ) - C j+2 (t) 
,Pij(x,t) := L ik (x)£ji(t)(x - Xi)(t-tj) - L i+2 {x)£ j+2 (t) 
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and 



Tl[u;xi + k,tj + i](x,t) := / (t - s)u tt (x l+k , s) ds 
Jtj+i 

+ (x-x i+ k) (t- s)u xtt (xi + k,s)ds 
Jtj+i 

(x- s)u xx (s,t) ds. 



The proof is very simple. In fact, (35 1 is a consequence of the following Taylor expansion for u 
with integral remainder: 



u(x,t) =u(x k ,ti) + (t-ti)u t (x k ,ti)+ / (t - s)u tt (x k ,s) ds 

Jti 

+ {x- x k ) \u x {x k ,t{) + (t- U)u xt {x k ,ti) + J^ (t - s)u x tt{x k ,s)ds 
+ J (x — s)u xx (s, t) ds 
and the identity ^ iji=< yLik{x)Cji(t) = 1. 

3.3 An estimate of \\p — Hhp\\p and some consequences 



Let us now prove that ( |31[ ) holds when the Ph are given by (33)-(34). 
Thus, let us fix p £ Pq and let us first check that 



p~ 2 \L{p- Il h (p))\ 2 dx dt -> as h = (Ax, At) (0,0). 
For each K k i £ Qh (simply denoted by K in the sequel), we write: 

JJ p- 2 \L( P -n h p)\ 2 dxdt<\\p- 2 \\ L ~ iK) JJ \L( P ~n hP )\ 2 dxdt 

< 3||p" 2 |U- W (^JJ^\(p-U h p) u \ 2 dxdt + JJ^\(a(x)(p-U hP ) x ) x \ 2 dxdt 

+ \\b\\ 2 L ^ (K) JJ K \ P -n hP \ 2 dxdty 

Using Lemma |3.3[ we have: 
JJ \p-U hP \ 2 dxdt 

= JJ y^X m ijPx(xj,tj) + nijPt{xi,tj) + PijPtx{%i,tj) +LiCjTl\p;Xi,tj]) 

K »J 

<l6\\Px\\l^ ( K)J2 jj \ m ij\ 2 dxdt + 16 \\Pt\\l-(K)J2 

i,3 K i,j 

+ 16|btx||z,°°qn ^ // \pij\ 2 dx dt + 16^2 \Li£jfc\p\%i,tj]\ 2 dxdt, 



(36) 



(37) 



dx dt 



\ 2 dxdt 



where we have omitted the indices k and I. 
Moreover, 



1,3 



|7e[p;x l+fc ,^ +i ]| 2 <\t- tj^Wpuixi, ■)\\ 2 L 2 ( t lttl+1 ) + \x- x^lt - t^Wp^Xi, -)\\ 2 L 2 (tutl+l) 

+ \X - Xi\ 3 \\p xx (- 1 t)\\ 2 L 2 (xk Xk+l y 
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Consequently, we get: 

E // \LiCjft\p;x i+ k,t j+ i]\ 2 dxdt 
ij JJk 

< sup \\p tt (x, Olll^t,,*.^) E // \U{x)Ci{t)\ 2 \t -tjfdxdt 
xe(x k ,x k+1 ) f . JJ/f 

+ sup HParttCrc, // iLiWCjittflt - tj\ 3 \x - xtf dx dt 

xE(x k ,x k + 1 ) JJK 

+ sup ||p xx (-,i)||| 2( Xk , Xk+l) J2 II \Li(x)£j(t)\ 2 \x-Xi\ 3 dxdt. 

After some tedious computations, one finds that 

^2jJ^\m ij \ 2 dxdt= ^1 (Aa .)3 At) J2jJ K \n> 3 \ 2 dxdt=^^Ax(Atf, 

W/ \ Pij \ 2 dxdt=-^-(Axf { Atf 



198450 



and 



E ff K \L i (x)CM 2 \t-t j \ 3 dxdt=^Ax(At)\ 

E |£i(*Ki(*)l 2 l*- *il 3 = ^)( Ax ) 4A< ' 
E |i 4 (x)£,(t)| 2 |x - ^| 2 |t - i,f dzdi = i ^ 5 (A a; ) 3 (Ai) 4 . 
This leads to the following estimate for any K = K^i £ Qh' 



JJ k \p H hP \ 2 dxdt <^{AxfAt\\ Px \\ 2 La 



(K) 



1664 A . A ,o n ll9 
+ 11025 Ax(A<) 

+ 9^ (Aa;)3(Ai)3||Pte|l ^W 

+ 1 ^Ax{Atf sup |b xx M)||| 2(xii , Xfe+l) 



We deduce that 



\p - Ll h p\ 2 dx dt <KiT\\p x f LOO{QT , (A 



KiT\\pt\\l-(Q T )W 



Qi 

, v T"-ll„ II 2 

(Qt) 

^ 2 T|M|^ (Qt) (Ax) 2 (At) 2 

^3|btt(-,t)|||2 (0 ,T;I,oo(0,l)) (Ax) 3 
^3 life {-Ml °°(0,T;L 2 (0,1) 

(At) 3 

^4||Px«(- ,t)||£2 ( o, T ;L~(0,l)) (A^) 2 (At) 2 
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for some positive constants K%. Hence, for any p £ Pq one has 

// \p-n hP \ 2 dxdt^0 as h->(0,0). 

JJqt 



Proceeding as above, we show that the other terms in (37 1 also converge to 0. Hence, (36) 
holds. 

On the other hand, a similar argument yields 

/ Po 2 a(l) 2 |(p- n h p) x \ 2 dxdt -> as h~^(0,0) 
Jo 



and, consequently, we find that (31 1 holds 



We can now use Proposition |3.1| and deduce convergence results for the approximate control 
and state variables: 



Proposition 3.2 Letph £ Ph be the unique solution to ( 30 1 7 where Py t is given by ( 33 ) — ( 34 1 . Let 
us set 

Vh ■= P~ 2 Lp h , v h := -pQ 2 a(x)p h . x \ x=1 . 

Then one has 

h ~ Vh\\mQ T ) -> and \\v - v h \\ L 2 {0 T} -> 0, 
where (y,v) is the solution to Q. □ 

3.4 A second approximated problem 

For simplicity, we will assume in this section that y\ £ C°([0, 1]). 

In order to take into account the numerical approximation of the weights and the data that we 
necessarily have to perform in practice, we will also consider a second approximated problem. It 
is the following: 

m h (p h ,q h ) = (£ h ,q h ) \/q h £ Phi Ph £ Ph, 
where the bilinear form TO/j(- , •) is given by 

m h (p h ,qh) ■= // n h (p~ 2 )Lp h Lq h dxdt + / a(l) 2 ir Ax (pQ 2 )p h q h dt 
JJqt Jo 

and the linear form £^ is given by 

(4,9):= / {^Axyo){x)q t {x,Q)dx- {■K Ax y u q{-,Q)) H - KH i 



(38) 



Here, for any function / £ C (Qt), Kh{f) denotes the piecewise linear function which coincides 
with / at all vertices of Qh- Similar (self-explanatory) meanings can be assigned to ttax{z) and 
n At (w) when z £ C°([0, 1]) and w £ C*°([0,T]), respectively. 

Since the weight p~ 2 is strictly positive and bounded in Qt (actually p~ 2 > 1), we easily see 
that the ratio ■Kh{p^ 2 ) I ' p~ 2 is bounded uniformly with respect to h (for \h\ small enough). The 



same holds for the vanishing weight 0jp(l, •) under the assumptions (28) 



As a consequence, it is not difficult to prove that (38) is well-posed. Moreover, we have 



Lemma 3.4 Let ph and ph be the solutions to (30) and (38), respectively. Then. 



\\Ph-Ph\\p <max 



1 



,7TAt(p ) 



1 



P "" Po 

+ Ci||7r A:E (y ) - y \\ L 2 + C 2 \\TTAx(yi) ~ Ui\\h-i, 

where G\ and C2 are positive constants independent of h. 



L°°(0,T) \\Ph\\P 



(39) 
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Proof: Since ph and ph respectively solve (30) and (38), one has: 



\\Ph -Ph\\p = m {ph -Ph,Ph - Ph) 

= m(ph,ph ~Ph) ~ rn h (Ph,Ph ~Ph) + (lh,Ph - Ph) - (l-,Ph - Ph) 

(p 2 -n h (p 2 ))Lp h L{p h -p h )dxdt+ I (p Q 2 -it^Pq ))a{l) 2 p x ,h{px,h-Px,h) dt 

JQ 

+ (lh)Ph - Ph) - (l,Ph - Ph) 

1 - ^j^) (p-'LphXp-'Liph - p h )) dx dt 



1 



n AtP 2 



a(l) (p Ph)p (Ph,x - Ph,x) dt 



+ / (TTAx(yo) - yo)(x) (Pt,h - pt,h)(x, 0) dx - (wAxiyi) - yi, {ph - Ph){x, 0)) H - KH i 
Jn 



In view of the definitions of the bilinear forms m{- , •) and m^(- , •), we easily find (39). 

Taking into account that ( [31] ) holds and 

^h{p- 2 



□ 



max 



P 



_ 2 ~ i \\L°°(Q T )i \\ n 1 



Po 



IU°°(0 : T)^ 



->o, 



we find that, as h goes to zero, the unique solution to (38), converges in P to the unique solution 



to (29): 



\\p-Ph\\p < \\p-Ph\\p + \\Ph ~Ph\\p -> o. 
An obvious consequence is the following: 



Proposition 3.3 Letph S Ph be the unique solution to (38), where Ph is given by (33)-(34). Let 
us set 

jjh ■= TTh{p~ 2 )Lp h , Vh ■= -^Ax{Po 2 )a{x)p htX \ x=1 . (40) 

Then one has 

II 2/ - yh\\v>(Q T ) -> and \\ v ~ W/i|U 2 (0,T) -> o, 
where {y,v) is the solution to rtSJ) . □ 



4 Numerical experiments 



We now present some numerical experiments concerning the solution of ( 38 ) , which can in fact be 



viewed as a linear system involving a banded sparse, definite positive, symmetric matrix of order 
AN x N t . We will denote by M. h this matrix. If {ph} stands for the corresponding vector solution 
of size 4N x N t , we may write {p h ,qh)p h = {M h {Ph}, {?/,}) for any q h € P h - 

We will use an exact integration method in order to compute the components of M.h and the 
(direct) Cholesky method with reordering to solve the linear system. 



After the computation of ph, the control Vh is given by (40). Observe that, in view of the 



definition of the space Ph, the derivative with respect to x of ph is a degree of freedom of {ph}', 
hence, the computation of Vh does not require any additional calculus. 

The corresponding controlled state jjh may be obtained by using the pointwise first equality 



(40) or, equivalently, by solving (24). However, in order to check the action of the control function 



Vh properly, we have computed jjh by solving with a C 1 finite element method in space and a 
standard centered scheme of second order in time. 
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Thus, let us introduce the finite dimensional spaces 

Zh = {z h eC 1 ([0,l]):z h \ [xi<Xi+Ax] €F 3<X Vi = l,...,N x } 

and Zoh = {zh G %h '■ z h{0) — %(1) = 0}. Then, a suitable approximation jjh of the controlled 
state y is defined in the following standard way: 

• At time t = 0, y ft is given by y h {-,0) = Pz h (yo), the projection of y on Z fe ; 

• At time t% = At, yh{',t\) G ^ft is given by the solution to 

" X (yh(x,ti) - yh(x,t Q ) - Atyi(x)) 



(Aty 



) (fa; 



+ / [a(x)yft, x (a;, i )</>z + 6(x, t )y h (x, t )4>] dx = 
Jo 

^ G Zha; Vh(0,ti) G Z/j, j)h(0,ti) = 0, = 



• At time t — t n — nAt, n — 2, • • • , A/j, i n ) solves the following linear problem: 

nl ($h(x,t n ) - 2y h (x,t n -i) + y h (-,t n - 2 ) 



(41) 



(Aty 



> dx 



+ / [a(x)y htX (x,t n -i)4> x + b(x,t n )yh(x,t n - 1 )4>] dx = 
Jo 

G Zfc ; Vh{0,t n ) G Z h) yh(0,t n ) = 0, y h (l,t n ) = v h (t n )- 



(42) 



This requires a preliminary projection of on a grid on (0, T) fine enough in order to fulfill 
the underlying CFL condition. To this end, we use the following interpolation formula: for any 
Ph G Ph and any 8 G [0, 1], we have: 



Ph^frtj + 6At) = [29 + 1)(0 - l) 2 p M (l, * 3 + At 0(1 - 0) 2 p Mt (l, t,-) 
+ # 2 (3 - 20)p M (l,t i+ i) + A^ 2 (0 - l)ph,xt(l,t i+ i) 



(43) 



for all t G [tj, ij+i]. 

We will consider a constant coefficient a(x) = clq = 1 and a constant potential b(x,t) = 1 in 



Q T - We will take T = 2.2, x = -1/20, /? = 0.99 and M = 1 - + /3T, so that (17) holds 



Finally, concerning the parameters A and s (which appear in (21 1), we will take A — 0.1 and s — 1 



Remark 6 Let us emphasize that our approach does not require in any way the discretization 
meshes to be uniform. □ 

4.1 Estimating the Carleman constant 

Before prescribing the initial data, let us check that the finite dimensional analog of the observabil- 



ity constant Co in (20) is uniformly bounded with respect to h when (17) is satisfied. We consider 



here the case a = 1 and 6=1. 



In the space Ph, the approximate version of (20) is 



(A h { Ph }, { Ph }) < C oh (M h {p h }, {ph}) V{ Ph } G P h , 
where Ah is the square matrix of order AN t N x defined by the identities 



(A h {ph}, Uh}) ■■= / (Ph,x(x, 0) qh,x(x, 0) + Ph,t{x, 0) qh,t(x, 0)) dx. 
Jo 
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Therefore, Coh is the solution of a generalized eigenvalue problem: 

C oh = max{ A : 3p h £ P h , p h ^ 0, such that A h {p h } = XM h {ph} }■ 



(44) 



We can easily solve ( 44 1 by the power iteration algorithm. Table [T] collects the values of C$h 
for various h = (Ax, At) for T = 2.2 and T = 1.5, with At — Ax. As expected, Coh is bounded in 
the first case only. The same results are obtained for At ^ Ax. 



Ax, At 


1/10 


1/20 


1/40 


1/80 


T = 2.2 


6.60 x 10~ 2 


7.61 x 10~ 2 


8.56 x 10~ 2 


9.05 x lO- 2 


T = 1.5 


0.565 


2.672 


17.02 


96.02 



Table 1: The constant Coh with respect to h. 

In agreement with Remark [4j we obtain the same behavior of the constant with respect to T 
for any s, in particular for s — leading to p = 1 and po = 1- 



4.2 Smooth initial data and constant speed of propagation 



We now solve ( 23 ) with a = 1 and smooth initial data. For simplicity, we also take a constant 
potential 6=1. 

For (yo,Ui) = (sin(7rx), 0), Table [2] collects relevant numerical values with respect to h = 
(Ax, At). We have taken At = Ax for simplicity but, in this finite element framework, any other 
choice is possible. In particular, we have reported the condition number n(M.h) of the matrix AAh, 
defined by 

K (X,) = |||M i ||| 2 |||^ 1 ||| 2 

(the norm |||Al/j|||2 stands for the largest singular value of M.h)- We observe that this number 
behaves polynomially with respect to h. 



Ax, At 


1/10 


1/20 




1/40 


1/80 




<M h ) 


3.06 x 10 8 


1.57 x 10 10 


6.10 x 10 11 


2.47 x 10 13 


\\Ph\\p 


1.541 x 10" 1 


1.548 x 10" 


-l 


1.550 x lO" 1 


1.550 x 10 


-l 


\\Ph-p\\p 


4.46 x 10~ 2 


1.45 x 10" 


2 


4.01 x 10~ 3 


8.38 x 10- 


4 


\\Vh L 2 (0,T) 


5.421 x 10" 1 


5.431 x 10 


-1 


5.434 x lO" 1 


5.434 x 10 


-3 


\\vh - v\\l2(0,T) 


2.39 x 10~ 2 


8.12 x 10" 


3 


2.48 x 10~ 3 


9.57 x 10- 


4 


\\Vh(- ,T)\\ LHOtl) 


1.80 x 10~ 2 


8.18 x 10" 


3 


1.64 x 10" 3 


5.85 x 10- 


4 


\\yt,h(- ,T)\\ H -i( 0A) 


3.06 x 10- 2 


8.25 x 10" 


3 


3.59 x 10~ 3 


1.93 x 10" 


3 


Table 2: (y (x), yi (x)) = 


(sin(7ra;), 0), 


a = 


1,6=1- T = 


2.2. 





Table [2] clearly exhibits the convergence of the variables pt and Vh as h goes to zero. Assuming 
that h — (1/160, 1/160) provides a reference solution, we have also reported in Table[2]the estimates 
1 1 jo — f>h\\p and \\v — Vh\\L 2 (o,T)- We observe then that 

\\p - p h \\ P = 0(h im ), \\v ~ v h \\ L 2 i0 .T) = 0(h 15e ). 

The corresponding state yh is computed from the main equation ([!]), as explained above, taking 
At = Ax/4. That is, we use Q with 6 = 0,1/4,1/2 and 3/4 on each interval [tj,tj+i\. We 
observe the following behavior with respect to h: 



\\yh(;T)\\ moA) = 0(h 171 ), \\y t , h (;T)\\ H -i {0 , T) = O^ 1 " 1 ), 
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which shows that the control Vh given by the second equality in ( 40 ) is a good approximation of a 
null control for ([I]). 

Figure [2]Teft displays the function p^ 6 P (the unique solution to (38)) for h — (1/80, 1/80). 



Figure [2]-Right displays the associated control ih- As a consequence of the introduction of the 
function 9s in the weight, we see that ih vanishes at times t = and t = T. Finally, Figure [3] 
displays the corresponding controlled state jjh- 




-1.6 -1.4 -1.2 

Iog10(h) 



Figure 1: log 10 \\p - ph\\p (*) and log 10 \\v - V h \\ L 2 {o,T) (°) vs. \og w {h). 



x 10" 3 




Figure 2: (y (x) , yi(x)) = (sin(7ra;), 0) and a = 1 - The solution p h over Q T (Left) and the 
corresponding variable Vh on (0,T) (Right) - h = (1/80, 1/80). 

Tableland Figures [4] and [i] provide the results for yo(x) = e _500 ( x_0 - 2 ) and yi(x) = 0. We still 
observe the convergence of the variables f>h, Vh and y^, with a lower rate. This is due in part to the 
shape of the initial condition y - Precisely, we get \\p— Ph\\p — ©(/i 1 ' 74 ), \\vh — u IIl 2 (o ,t) = O(h om ), 
\\yhi;T)\\mo,i) = 0(h 135 ) and \\yt.h{;T)\\ H -^ .T) = 0{h l u ). 
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o o 



Figure 3: [yo{x),y\{x)) = (sin(7ra;), 0) and a = 1 - The solution jjh over Qt - h = (1/80,1/80). 



Ax, At 


1/10 




1/20 




1/40 




1/80 




1/160 




\\Ph\\p 


4.38 x 10~ 


-2 


3.95 x 10~ 


-2 


4.20 x 10- 


-2 


4.31 x 10- 


-2 


4.33 x 10- 


-2 


WPh ~P\\p 


1.80 x 10" 


-1 


6.30 x 10- 


-2 


1.66 x 10- 


-2 


2.78 x 10- 


-3 






\\Vhh*(0,T) 


1.48 x 10~ 


-1 


1.33 x 10- 


-1 


1.53 x 10- 


-1 


1.64 x 10- 


-1 


1.67 x 10- 


-1 


\\vh - v\\li{q,t) 


9.81 x 10" 


-2 


6.28 x 10" 


-2 


3.80 x 10- 


-2 


1.11 x 10- 


-2 






\\Vh{- ;7 1 )|| i 2 (0 1} 


1.09 x 10~ 


-1 


7.67 x 10- 


-2 


3.70 x 10- 


-2 


1.11 x 10- 


-2 


1.87 x 10- 


-3 


\\Vt.h{- > T )l!ff- 1 (0,l) 


1.36 x 10~ 


-1 


8.82 x 10- 


-2 


5.16 x 10- 


-2 


1.76 x 10- 


-2 


2.82 x 10- 


-3 



Table 3: (yo(x), yi (x)) = (p-500^-0.2) 2 > y an d a = 1,6= 1 - T = 2.2. 



x 10" 4 




Figure 4: (y (x) , yi(x)) = (e 500 ^ 2 ' 2 ,0) and a = 1 - The solution f>h over Qt (Left) and the 
corresponding variable Vh on (0,T) (Right) - h = (1/80, 1/80). 
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Figure 5: (yo(x),y 1 (x)) = ( e - 500 ( a; - - 2 ) 2 , 0) and a = 1 -The solution y h over Q T - h = (1/80, 1/80). 



4.3 Initial data (yo,yi) G H 1 ^, 1) x L 2 (0, 1) and constant speed of propa- 
gation 

Let us enhance that our approach, in agreement with the theoretical results, also provides conver- 
gent results for irregular initial data. We take a continuous but not differentiable initial state yo 
and a piecewise constant initial speed y\. 

y (x) = x l[ ,i/2]O) + (1 - x) l]i/ 2 ,i]( a: ) I 2/1 fa) = 10 X 1 [1/5,1/2] (x). (45) 

The other data arc unchanged, except b, that is taken equal to zero. 

Observe that these functions remain compatible with the C 1 finite element used to approxi- 
mate p, since y$ and yi only appear in the right hand side of the variational formulation and ttaxVo 

18). The unique difference is that, once ph and ih are known, y^ must 
using a C° (and not C 1 ) spatial finite element method. 
Recall however that these initial data typically generate pathological numerical behavior when 
the usual dual approach, based on the minimization of (J©]) , is used. 

Some numerical results are given in Table [4] and Figures [6] and [7] As before, we observe the 
convergence of the variable ph and therefore Vh and y/,asft-> 0. We see that \\ph — p\\p = 0(h 1AS ) 
and \\vh — ^||l 2 (o,i) = 0(h 123 ). In particular, we do not observe oscillations for the control or the 
functions ph and p>h,t at the initial time. 



and ttaxVi make sense; see ( 
be computed from (41)-(W2 



Ax, At 


1/10 




1/20 




1/40 




1/80 




1/160 


\\Ph\\p 


3.16 x 10" 


1 


2.89 x 10" 


2 


2.73 x 10" 


2 


2.65 x 10" 


2 


2.61 x 10- 1 


\\Ph-p\\p 


1.12 x 10" 


1 


4.62 x 10" 


2 


1.70 x 10" 


2 


5.12 x 10" 


-3 




H*/ilU*(0,T) 


1.23 




1.11 




1.05 




1.02 




1.004 


\\i>h - «IU 2 (0,T) 


2.52 x 10" 


1 


1.25 x 10" 


1 


5.57 x 10- 


2 


1.90 x 10" 


2 




\\y~h{- ,T)\\ L 2 {0A) 


1.09 x 10- 


1 


5.40 x 10" 


2 


2.20 x 10" 


2 


1.09 x 10" 


-2 


6.20 x 10~ 3 


\\yt,h(- : T )llff- 1 (0,l) 


7.25 x 10" 


2 


4.62 x 10" 


2 


2.85 x 10" 


2 


5.12 x 10" 


-3 


6.75 x 10~ 3 



Table 4: (y , yi) given by (45 1 and a = 1 - T = 2.2. 



4 NUMERICAL EXPERIMENTS 



21 




Figure 7: (2/0,2/1) given by (45 1 and a = 1 -The solution y/i over - h = (1/80, 1/80). 
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4.4 Discontinuous initial data y and constant speed of propagation 



The method also provides convergent results for data yo only in L 2 (0, 1). We consider the following 
initial condition: 

Vo{x) = l[ .5,o.7](x), yi(x) = Q. (46) 

The other data are unchanged with respect to Section [4~3] This leads to pathological numerical 
behavior when other frequently used dual methods are employed (we refer to |24j ) . Some numerical 
results are given in Table [5] and Figure [8j Once again, the convergence of the variable ph and 
therefore Vh and ijh as h — > is observed. 



Ax, At 


1/10 


1/20 


1/40 


1/80 


1/160 


\\Ph\\p 


1.01 x lO" 1 


1.00 x lO" 1 


9.71 x 10^2 


9.53 x lO- 2 


9.47 x 10~ 2 


II«/i||l 2 (0,T) 


3.42 x 10" 1 


3.27 x 10" 1 


3.19 x 10" 1 


3.14 x 10" 1 


3.14 x 10" 1 


\\Vh(- ■T)\\l^(o,i) 


1.24 x lO" 1 


9.27 x 10~ 2 


7.26 x 10~ 2 


5.88 x 10~ 2 


3.12 x 10" 2 


\\yt,h(- ,T)\\ H -i( 0t i) 


1.55 x 10" 1 


1.16 x lO" 1 


1.06 x lO" 1 


7.13 x 10" 2 


6.02 x 10~ 2 



Table 5: (y ,yi) given by (46) and a = 1 - T = 2.2 




4.5 Non constant smooth speed of propagation 

Finally, let us consider a non-constant function a — a(x) (we refer to [5D] for the dual approach 
in this case). In order to illustrate the robustness of our method, we will take a coefficient a € 
(^([O.l]) with 

r 1 x G [0, 0.45] 

a(x) = I €[1.,5.] {a'(x)>0), x € (0.45, 0.55) (47) 
[5 x e [0.55, 1] 

so that condition (17) is equivalent to T > 2(1 + l/20)-\/5 ~ 4.69 (taking again xq = —1/20). In 
order to reduce the computational cost, we take as before T = 2.2 and we still observe that the 
constant Coh in (44) is uniformly bounded. 
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We take again (y (x), yi(x)) = ( e - 50 °( x - a - 2 ) 2 : o) and 6 = 0. Table [fj] illustrates the convergence 
of the approximations with respect to h. Figures [9] and 10 depict for h = (1/80, 1/80) the functions 
Ph, Vh and jjh- In particular, in Figure[lOj we can observe the diffraction of the wave when crossing 
the transitional zone (0.45,0.55). 




o 



Figure 9: (y (x),yi(x)) = (e 500 ( x ° ,0) and a given by (47) - The solution p h over Q T (Left) 



and the corresponding variable on (0,T) (Right) - h = (1/80, 1/80). 



Ax, At 


1/10 




1/20 




1/40 




1/80 




1/160 




\\Ph\\p h 


3.87 x 10" 


2 


3.44 x 10" 


2 


3.75 x 10- 


2 


3.85 x 10- 


2 


3.86 x 10" 


2 


\\Ph-p\\p h 


1.25 x 10" 


1 


5.75 x 10- 


2 


2.64 x 10- 


2 


1.01 x 10- 


-2 






IMU»(0,T) 


7.74 x 10- 


2 


6.53 x 10- 


2 


9.16 x 10- 


2 


1.01 x 10- 


1 


1.03 x 10- 


-1 


114 - "IU 2 (0,T) 


5.07 x 10- 


1 


4.17 x 10- 


2 


2.03 x 10- 


2 


4.86 x 10- 


-3 






l|yh.(-.^)llL 2 (o,i) 


1.09 x 10- 


1 


7.89 x 10- 


2 


1.81 x 10- 


2 


1.16 x 10- 


-2 


1.71 x 10" 


3 


\\yt,h(- ,T)\\ H -i(o,i) 


1.01 x 10- 


1 


8.39 x 10- 


2 


4.81 x 10- 


2 


7.52 x 10- 


-3 


1.55 x 10- 


-3 



Table 6: (y (x),yi(x)) = ( e -soo(x-o.2) ,o) and a given by (47| _ T = 2 2 



5 Further comments and concluding remarks 

Let us begin this section with some general considerations on the use of Carleman weights that 
serve to justify our approach: 

(i) The search of a control minimizing J in ([8]), where y is involved, is very appropriate from 
the numerical viewpoint. As shown in Section [2j the explicit occurrence of the state variable 
y leads to an elliptic problem in Qt, that is easy to analyze and solve (at this level, the 
particular choice of the weight is less important). This approach does not require the dis- 
cretization of the wave operator, as for usual dual approachs; therefore, it does not generate 
any spurious oscillations and leads to numerical well-posedness. This is an important feature 
of the approach. 



(ii) The Carleman weights provide regularity of the solution to ( 23 1 and therefore allows to 
derive estimates of the errors \\p — Ph\\p in term of h = (Ax, At). This will be detailed in a 
forthcoming work. 
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Figure 10: y (x) = e -soo(x-o.2) and a giyen by (47) _ The solution yh over Q T -h = (1/80, 1/80) 



(iii) The process can be viewed as a first step for the numerical controllability of semi-linear 
problems: if we just apply a fixed-point argument, we will find at each iterate a linear 
equation with non-regular coefficients depending on x and t for which the present approach 
is adequate. 

(iv) In our numerical experiments we have not found any essential difference for small or large s 
or A: this is in full agreement with Remark [4] and Section |4~Tj 

(v) Furthermore, we mention that the approach has been considered by analogy with the similar 
analysis in references [T2] to [TS], dealing with heat equations. 



5.1 Primal versus dual approach (I): analogies 



The solution to the variational formulation (|23j) is also the unique minimizer of the functional I, 
with 

I(P)-=1 [[ p- 2 \L P \ 2 dxdt+ l - [ pfa{lf\p x {l,-)\ 2 dt 

1 Jo (48) 



J(-0) 



(h 



y o (x)p t (x,0) dx + (yi,p{; Q)) H -i, H i. 

This is similar to the conjugate functional J* in Actually, we notice that J*(/i, 4>o, <t> 
for all (/x, 4>q,4>i) £ L 2 (Q T ) X H. 

Therefore, the extremal problems ^ and ( 48 1 are connected to each other having ^ as starting 
point. The problem (23), deduced from the primal approach belongs to the framework of elliptic 



variational problems in two dimensions and is well tailored for a resolution with finite elements. The 
dual problem ([9]) is of hyperbolic nature: the time variable is kept explicitly and time integration 
is required. 

Note that we may also derive the optimality conditions for J* (as we did in Section [2] for J): 
this leads, at least formally, to the problem (23). 



We also mention [26] where a (different) variational approach is introduced. 
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5.2 Primal versus dual approaches (II): discrete properties 



The variational approach used here leads to satisfactory convergence results, in particular the 
strong convergence of the approximate controls Vh towards a null control of the wave equation. 
This relies in a fundamental way on the fact that we work in a subspace Ph of P. Indeed, this 



allows to write directly the Carleman estimate in Ph and get that the function / (given by (48)) is 
uniformly coercive with respect to the discretization parameter h. 

On the other hand, notice that no wave equation has to be solved in order to compute the 
approximations Vh- For each h, once Vh is known, we must solve the wave equation, in a post- 
treatment process, to compute the corresponding state jjh (recall that, actually, this may be avoided 
by using directly the optimality condition y = —p~ 2 Lp). 

This is in contrast with the dual approach. Indeed, the minimization of J* by an iterative pro- 
cess requires the resolution of wave equations, through a decoupled space and time discretization. 
As recalled in the introduction, this may lea to numerical pathologies (the occurrence of spurious 
high frequency solutions) and, therefore, needs some specific numerical approximations and tech- 
niques. We mention the work [5], where the authors prove, in a close context and within a dual 
approach, a weaker uniform semi-discrete Carleman estimate with an additional term in the right 
hand side, necessary to absorb these possibly spurious high frequencies (see [5], Theorem 2.3). 

Notice that the computed Vh are not a priori null controls for discrete systems (associated to 
the wave equation ([!])), but simply approximations of the control v furnished by the solution to ([8]). 
If one wants to go further in the comparison, it can be said that the primal approach aims to first 
compute the control for ([!]) and then approximate it, while the dual classical method aims first to 
discretize and then control the corresponding finite dimensional system. 

Let us also observe that the (primal) approach in this paper is relatively easy to implement. 
In practice, the resolution is reduced to solve a linear system, with a banded sparse, symmetric 
and definite positive matrix, for which efficient direct LU type solvers are known and available. 
Furthermore, we may want to adapt (and refine locally) the mesh of Qt in order to improve 
convergence and such adaptation is much simpler than in the dual approach, where t is "conserved" 
as a time variable. For additional considerations, see also [5] and [5]. 

5.3 Mixed formulation and C°-approximation 

The approach can be extended to the higher dimensional case of the wave equation in a bounded 
set ft C H. , with N > 2. However, the use of C 1 -finite element is a bit more involved. Arguing 
as in " 



j, we may avoid this difficulty by introducing a mixed formulation equivalent to (23). 



The idea is to keep explicit the variable y in the formulation and to introduce a Lagrange 
multiplier, associated to the constraint p 2 y + Lp = (see (22)). We obtain the following mixed 
formulation: find (y,p, A) € Z X P X Z such that 



where 



p 2 yydxdt + [ p 2 a(l) 2 Px {l,t)p x (l,t)dt + [[ \{p 2 y + Lp) 

JO JJQt 



y {x)p t (x,0)dx - {yi,p{-,ti)) H -\m V(y,p) e Z x P, 



X(p 2 y + Lp) dxdt = Q VA e Z, 



Z = L 2 (p 2 ;Q T ) := {z e L\ oc {Q T ) : [[ p 2 \z\ 2 dx dt < +oo }. 

JJQt 



dx dt 



(49) 



Taking advantage of the global estimate (16), we may show, through an appropriate inf-sup condi- 



tion, that ( 49 ) is well-posed in Z x Px Z. Moreover, the approximation of this formulation may be 
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addressed using C°-finite element, which is very convenient. The approximation is non-conformal. 
More precisely, the variable p is now sought in a space Rh of C°-functions that is not included in 
P. 



At the discrete level, (49 1 reduces the controllability problem to the inversion of a square, 
banded and symmetric matrix. Moreover, as before, no wave equation has to be solved, whence the 
numerical pathology described above is not expected. However, since the underlying approximation 
is not conformal (this is the price to pay to avoid C 1 finite elements), a careful (and a priori not 
straightforward) choice for R^ has to be done in order to guarantee a uniform discrete inf-sup 
condition. The analysis of this point, as well as the use of stabilized finite elements, will be 
detailed in a future work. 



5.4 Extensions 

The approach presented here can be extended and adapted to other equations and systems. What 
is needed is, essentially, an appropriate Carleman estimate. 

In particular, we can adapt the previous ideas and results to the inner controllability case, 
i.e. the null controlability of the wave equation with distributed controls acting on a (small) sub- 
domain oj of (0, 1). Furthermore, using finite element tools, we can also get results in the case 
where the sub-domain w varies in time, that is non-cylindrical control domains qx of the form 

q T = { (x,t) e Q T : gi(t) <x<g 2 (t), t £ (0, T) }, 

where g± and gi are smooths functions on [0, T], with < gx < gi < 1. This opens the possibility 
to optimize numerically the domain qj>, as was done in a cylindrical situations in I25j (see also 

mi 

Let us finally mention that many non-linear situations can be considered through a suitable 
linearization and iterative process. We refer to [141 I15j for some ideas in a similar parabolic 
situation. 



A Appendix: On the proof of Theorem 2.1 

We first prove a global Carleman estimate for functions w satisfying vanishing initial and final 
conditions. In what follows, L stands for the operator given in ([2) with 6=0. It is easy to 
check that, if the estimate ( 16 ) holds in this particular case, then the same estimate holds for any 
potential b £ L°°((0, 1) x (— T, T)). 

Theorem A.l With the notation of Section^ let xq < be a fixed point, let <f> and ip be the 
weight functions defined by ( 12 I -( 14 1 and let a £ A(xq, Oq) with ao > 0. Then there exist positive 
constants sq and M, only depending on Xq, oq, ||a||c 3 ([o.i]) an d T such that, for all s > sq, one 
h as: t l T 1 

s[ f e 2sv (\v t \ 2 + \v x \ 2 ) dxdt + s 3 [ f e 2sip \v\ 2 dxdt 

pT pi pT y° u ) 

<M / e 2sip \Lv\ 2 dxdt + Ms e 2sv \v x (l, t)\ 2 dt 



J-TJO J-T 

for any v G L 2 (—T, T; Hq(0, 1)) satisfying Lv £ L 2 ((Q,1) x (— T, T)), v x (l, •) £ L 2 (-T,T) anc 



v(-,±T)=v t (-,±T)=Q. 
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The proof of this result follows step-by-step the proof of Theorem 2.1 in [3j. However, since the 
argument provides conditions on the set of admissible a and, to our knowledge, these conditions 
have not been stated in this form before, we provide here the detailed proof. 

Proof: Let us introduce w — e stp v and let us set 

Pw := e s <?L(e-^w) = e s <? ((e-^w) tt - (a( e -^w) x ) x ) . 
After some computations, we find that Pw = P\w + P 2 w + Rw, with 
Piw = w tt - (aw x ) x + s 2 X 2 ip 2 w (\4> t \ 2 - a\ip x \ 2 ) 

P 2 w = (a - l)s\(pw (iptt - {a^x)x) - s\ 2 ipw (|?/>t| 2 - a# x | 2 ) - 2sXtp (ip t w t - aip x w x ) 
R w = -asXcpw (ipu - (aip x ) x ) , 

where the parameter a will be chosen below. 
Recall that 

i){x,t)= \x - x \ 2 - Pt 2 + M , tp(x,t) = e^ x > t} 

and 

#M)>1 V(.x,<) g (0,1) x (-T,T). 

In this proof, we will denote by M a generic positive constant that can depend on xq, ao, ||a||c 3 ([o,i]) 
and T. 

As in the constant case a = 1, the first part of the proof is devoted to estimate from below the 
integral 



T 3 

1 = / f (Piw)[P 2 w)dxdt= V hj. (51) 
J-tJci id=1 



By integrating by parts in time and/or space, we can compute the integrals Iij in (51). We 
obtain: 

T i-l 

w tt V>w(ip tt - (atp x ) x ) dx dt 



hi = {a- l)sX f f 

J-T JO 



(l-a)sX / Lp\w t \ 2 (^ tt - (aip x ) x ) dxdt 

J-T JO 

(1 — a) f T f 1 

— sX 2 / <p\w\ 2 i/) tt (ip tt - (aijj x ) x )dxdt 

z J-T Jo 

-i^lsA 3 f [\\w\ 2 \ijt\ 2 ^u-(a^ x ) x )dxdt, 

1 J-T Jo 



1x2 = -sX 2 I I wtt <pw(\ip t \ 2 - a\ip x \ 2 ) dxdt 

J-T Jo 

= sX 2 f f ip\w t \ 2 {\iPt\ 2 -a\i/j x \ 2 )dxdt-sX 2 [ f V \w\ 2 \M 2 dxdt 

J-T Jo J-T Jo 

3sA 3 f T f 1 sX 3 f T f 1 

/ / ip\w\ 2 \ipt\ 2 i>ttdxdt + —- / / (p\w\ 2 a\ip x \ 2 iptt dxdt 
1-tJo z J-tJo 



2 

sX 4 r T r 1 

2 J-T JO 



ip\w\ 2 \^ t \ 2 (\H 2 -a\^x\ 2 )dxdt 
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and 



Ji3 = -2sX / Wtt<p(iptWt-aipxWx))dxdt 
J-t Jo 

= sX / if\w t \ 2 ipttdxdt + s\ 2 / / Lp\w t \ 2 \ip t \ 2 dxdt 

J-T Jo J-T Jo 

+sX / / (p\w t \ 2 {aijj x ) x dx dt + sX 2 / / ip\w t \ 2 a\tp x \ 2 dx dt 

J-T Jo J-T Jo 

— 2sX 2 / / (paipxiptWxWtdxdt. 



Also, 



'22 



and 



t Jo 



I21 = {l-a)sX / (aw x ) x ipw{iptt - {aip x ) x )dxdt 

J-T Jo 

= -(l-a)sX / tpa\wx\ 2 (iptt- (aipx)x)dxdt 

J-T Jo 



-T JO 
r T /■! 



(1 — ct) f f 

H = — sA2 / / ^M 2 (a^)x(^tt - (aip x ) x ) dx dt 

1 J-t Jo 

(1 — a) f T f 1 

+ 2 sX 3 J (pa\w\ 2 \^ x \ 2 (i;tt - (aip x ) x )dxdt 



-T JO 
fT r l 

\2 



(l-a)sA / / <p\w\ aip x (aip x ) xx dxdt 
t Jo 

f-T t l 



r T r l 

\2 



(1 — Ot) I f 

sX / ip\w\ 2 {a x (aip x ) xx + a(aip x ) xxx ) dxdt, 

J-T Jo 

sX 2 I I {aw x )x^w(\tp t \ 2 - a\ip x \ 2 ))dxdt 
J-t Jo 

-sX 2 / ipa\w x \ 2 (\ip t \ 2 - a\ip x \ 2 ))dxdt 
J-t Jo 

sX 2 f T f 1 

-— j J tp\w\ 2 ((\a x \ 2 + aa xx )\^x\ 2 + Aaa x ipxipxx + 2a(aip x ) x ip xx ) dxdt 
ip\w\ 2 (atp x )x{\ipt\ 2 ~ a\ip x \ 2 )) dxdt 



SA 

~2 



T JO 



sX 4 f T f 1 

+ / ¥\w\ 2 a\^x\ 2 (\*Pt\ 2 -a\4>x\ 2 ))dxdt 

1 J-T Jo 

-sX 3 / tp\w\ 2 aip x (a x \ipx\ 2 + 2aip x ip xx ) dxdt 
J-t Jn 



2sX / (aw x ) x tp(ip t w t - aipxWx) dxdt 

J-T Jo 

sX / ipa\w x \ 2 (ip tt + aip xx )dxdt 

J-T Jo 



I23 — 



+sX 2 / / ipa\w x \ 2 (iV'tl 2 + a l^a;| 2 ) dxdt — 2sX 2 / / (paipx-iptWxWtdxdt 

J-T Jo J-T Jo 

-sX J [a(l) 2 | Wx (l, t)\ t)^ x {l, t) - a(0) 2 \w x (0, t)\ V(0, t)i/> x (0, t)] dt. 
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Finally, 



J 3 i = (a ~ l)s 3 A 3 f [ <p 3 \w\ 2 {\i> t \ 2 -a\i> x \ 2 ){i> tt -{a^ x ) x )dxd^ 

J-T Jo 

/ 32 = -s 3 A 4 f {\ 3 \w\ 2 {\^ t \ 2 - a\^ x \ 2 ) 2 dxdt 

J-T Jo 



and 



hs = -2s 3 A 3 f [ v 3 w(\ip t \ 2 -a\ip x \ 2 )(ij t w t -aip x w x )dxdt 

J-T Jo 

= s 3 X 3 / / ^ 3 \w\ 2 (\ij t \ 2 -a\^ x \ 2 )(iP u -(aiP x ) x )dxdt 

J-T Jo 

+2s 3 X 3 [ f cp 3 \w\ 2 (\i>t\ 2 ip tt + aa x ip x \ip x \ 2 + a 2 \ip x \ 2 ip xx ) dxdt 

J-T Jo 

+3s 3 A 4 / / if 3 \w\ 2 (\ij: t \ 2 -a\iP x \ 2 ) 2 dxdt. 

J-T Jo 

Gathering together all terms for i, j € {1,2,3}, we obtain 

pT r \ 

1 = 



f [ (Piw) (P 2 w) dx dt 

J-T Jo 



= sX / ip\w t \ 2 (2ip tt - a(ip u - (atp x ) x )) dxdt 

J-T Jo 

+ sX / tpa\w x \ 2 (a(iptt - (aip x ) x ) + 2(aip x ) x - a x ip x ) dx dt 

J-T Jo 

+ 2sX 2 / ip (\w t \ 2 \ip t \ 2 - 2aip x ip t w x w t + a 2 \w x \ 2 \ip x \ 2 ) dxdt 
J-t Jo 

+ 2s 3 X i f [\ 3 \w\ 2 (]ip t \ 2 - a\ip x \ 2 f dxdt 

J-T Jo 

+ s 3 X 3 [ [ ip 3 \w\ 2 (\2ip t \ 2 tptt + aa x ip x \ip x \ 2 + 2a 2 \i> x \ 2 t/j xx ) dxdt 

J-T Jo 

+ as 3 X 3 [ [ cp 3 \w\ 2 (\ip t \ 2 -a\ip x \ 2 )(iJu-(aip x ) x )dxdt 

J-T Jo 

-sX [ (o(l) 2 |u; a .(l,«)|V(l.*)V'*(l,t)-o(0) a |i« a! (0,t)| 2 ^(0,t)V» a ,(0 ) t)) dxdt 

J-T 

+ x , 

where Xq is the sum of all "lower order terms" : 



\X \ < MsX 4 [ [ 

J-T Jo 



T pi 

3|„„|2 



<p \w\ dx dt. 



Let us analyze the high order terms arising in the previous expression of I. First, remark that 
sX 2 f [ ip(\w t \ 2 \A\ 2 -2a^ t w x w t + a 2 \w x \ 2 \^ x \ 2 ) dxdt>0. (52) 

J-T Jo 



Secondly, notice that, under the assumption a € A(xQ 1 a ), if (3 satisfies (11), we can choose a 



in such a way that the terms of order sX are positive. Indeed, we have in this case 
- a(x) - (x - x )a x (x) < f3 < a(x) + -(x - x )a x (x) Vx E [0, 1], 
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whence 

2/3 < 2a(x) + (x - x Q )a x (x) ^ g 1 

/3 + a(x) + (x — xo)a x (x) j3 + a(x) + (x — xq)cl x (x) 



Let a satisfy 



2/3 \ / 2<z(x) + (x - x )a x (x) 

sup r — -— < a < vox 1 



[0)1 ] \/3 + a(x) + (x - x Q )a x {x)J [o,i] V- 5 + a ( :r ) + ( x ~ xo)a x {x) _ 

Then, an explicit computation of the derivatives of V> shows that 

2ip tt - a(t/j tt - (aip x ) x ) > and a(V>tt - (aip x ) x ) + 2(aip x ) x - a x ip x > in [0, 1] x [-T, T] 
and, consequently, 

sX / / ip|u> t | 2 (2t/j tt - a(if) tt - (aip x ) x )) dxdt 

J-T JO 

+ sX / tpa\w x | 2 (a(ip tt - (a%l> x ) x ) + 2(aip x ) x - a x tp x ) dxdt 
J-t Jo 

>Msx[ [ <p\w t \ 2 dxdt + MsX [ f ip\w x \ 2 dxdt. 

J-T JO J-T Jo 

The remaining terms in / can be written in the form 

2s 3 A 4 f (\ 3 \w\ 2 {\^ t \ 2 - a\^ x \ 2 ) 2 dxdt 
J-t Jo 

+ s 3 X 3 ( [ ip 3 \w\ 2 (2\ip t \ 2 if>tt + aa x i) x \ip x \ 2 + 2a 2 \i/j x \ 2 iP xx ) dx dt 

J-T Jo 

+ as 3 X 3 [ [ <p 3 \w\ 2 (\tp t \ 2 -a\t/j x \ 2 )(i) tt -(aip x ) x )dxdt 

J-T Jo 

3 A 3 / / ip 3 \w\ 2 F x (x,Y(x,t))dxdt 7 



— s 

J-T JO 

where Y := \ipt\ 2 — a|'0£c| 2 an d 

( F x (x, Y) := 2AF 2 + (2ty« + - (a{x)^ x ) x )) Y 

+ a(x)\4> x \ 2 (2ip tt + a x (x)ip x + 2a(x)ip xx ) 
= 2XY 2 + (4/3 + a(2/3 + a(x) + (x - x )a x (x))) Y 
+ 8a(x)(x — xo) 2 (— 2(3 + 2a(x) + (x — xo)a x (x)). 

Since F\ is polynomial of the second degree in Y , one has 

F\(x, Y) > 8a(x)(x - x ) 2 (-2/3 + 2a(x) + (x - x )a x (x)) 

- —[4/3 + a(20 + a(x) + (x - x )a x (x))] 2 



for all x € [0,1] and Y € R. Therefore, if /3 satisfies |TT| ), for A large enough (depending on x 
and ||a||c 3 ([o.i]))j we obtain: 

s 3 \ 3 [ [ Lp 3 \w\ 2 F x {X)dxdt> Ms 3 X 3 [ [ tp 3 \w\ 2 dxdt. (53) 

J-T Jo J-T Jo 



Putting together the estimates (|52[)-(53), the following is found: 

/ [ (Piw) (P 2 w) dx dt > MsX [ I ip (\w t \ 2 + \w x \ 2 ) dxdt 

J-T Jo J-T Jo 



+Ms 3 X 3 [ [ tp 3 \w\ 2 dxdt (54) 

J-T Jo 

fT pi 

-MsX I \w x (l, t)\ 2 dx dt — MsA 4 

T 



[ [ ip 3 \w\ 2 dxdt. 

J-T Jo 
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On the other hand, recalling the definition of P, Pi, P 2 and R, we observe that 



/ / {\P lW \ 2 + \P 2 w\ 2 ) dx dt + 2 / f (P lW )(P 2 w)dxdt 
It is not difficult to sec that there a exists M such that 



T P l 



T r l 



\Pw-RwV dxdt 



T JO 



T r l 



\Pw - Rw\ 2 dx dt < M 



l-T JO 

In particular, we have 

fT pi 



T pi 



\Pw\ z dxdt + Ms 2 X 2 



T JO 



T pi 



f 2 \w\ 2 dxdt. 



/i pi pT pi pi pi 

/ (P 1 w)(P 2 w) dx dt < M / / \Pw\ 2 dxdt + Ms 2 X 2 / Lp 2 \w\ 2 dxdt 
-T JO J-T JO J-T JO 



-T JO 



T pi 



(55) 



and combining ( 54 ) and ( 54 ) we obtain: 



sA 



T pi 
-T 



f <f (\w t \ 2 + \w x \ 2 ) dx dt + s 3 X 3 [ I (p 3 \w\ 2 dxdt 

J-T JO 

M [ [ \Pw\ 2 dxdt + MsX [ \w x (l,t)\ 2 dt 

J-T Jo J-T 



MsX 4 



T pi 



ip 3 \w\ 2 dxdt + Ms 2 X z 



T JO 



T pi 



(p 2 |w| 2 dx dt. 



T JO 



Obviously, the last two terms in the right hand side can be absorbed by the second term in 
the left for s large enough. Therefore, there exists sq > 0, only depending on xq, o,q, ||a||c 3 ([o,il) 
and T, such that, for all s > sq, one has: 



T pi 



X I I ip (\w t \ 2 + \w x \ 2 ) dxdt + s 3 X 3 



T JO 



< M 



T pi 



ip 3 \w\ 2 dx dt 



T JO 



T pi 



T JO 



\Pw\ 2 dx dt + MsX I \w x {l,t)\ 2 dt. 

-T 



(56) 



Since w = ve sv and Pw = e stp Lv, we can easily rewrite (56) in the form (50 1. 
This ends the proof. 



□ 



In the remaining part of the Appendix, we will use the Carleman estimate ( 50 ) to prove The- 
orem 

Thus, let us assume that ^ holds , w E L 2 (-T, T; H%(0, 1)), Lw € L 2 ((0, 1) x (-T, T)) and 
w x (l, •) £ L 2 (-T,T). Thanks to §L5§, there exists r] € (0,T) and e > 1 such that 

(1 - e)(T - rj)P > max a(x) 1/2 (x - x ). 

[0,1] 

Moreover, simple computations show that, for every t £ (-T, — T + rj) U (T — rj, T), the func- 
tion , t) satisfies: 



(1 — e) min \ib t (x, t)\ > max aixV^libxix, t)\ 

[0,1] [0,1] 

max ib(x. t) < min ib(x, 0). 

. [0,1] [0,1] 

Let x £ (IK) a cut-off function such that < x < 1 and 

m fl, i{\t\<T- V 



(57) 
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Then we can apply Theorem A.l to the function w :— x w : whence the following Carleman estimate 
holds t 1 t 1 

si I e 2sip (\w t \ 2 + \w x \ 2 ) dx dt + s 3 [ [ e 2sip \w\ 2 dxdt 

<M f [ e 2sv \Lw\ 2 dxdt + Ms f e 2sv \w x (l, t)\ 2 dt. 

J-TJO J-T 



Since Lw = x^ w + Xttw + 2\tWt, we deduce from (57) that 



r T-n pi pT-n pi 

s / e 2sv (\w t \ 2 + \w x \ 2 ) dxdt + s 3 / / e 2s *\w\ 2 dxdt 

J-T+n JO J-T+n Jo 

<M f ( e 2sv \Lw\ 2 dxdt + Ms [ e 2stp \w x (l, t)\ 2 dt (58) 

J-T Jo J-T 

p-T+n pi pT pi 

+M / e 2sv (\w t \ 2 + \w\ 2 )dxdt + M / / e 2sv (\w t \ 2 + \w\ 2 ) dx dt. 

J-T Jo Jt-ti Jo 

Let us denote by E s = E s (t) the energy associated to the operator L, that is, 

E.(t) := \ J o e 2s * {\w t \ 2 + a\w x \ 2 ) dx. (59) 
Then, the argument employed in the proof of Theorem 2.5 in [2] (using the modified energy given 



by (59 1) can be used to deduce (161 from (58 1 
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